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Summary 

Guidelines are presented for use of the computer program PANCOR to assess the interference due\ 
to tunnel walls and model support in a slotted wind tunnel test section at subsonic speeds. Input data \ 
requirements are described in detail and program output and general program usage are described. The 
program is written for effective automatic veclorization on a CDC CYBER 200 class vector processing 
system. 

Introduc t ion 

PANCOR is a computer program for assessing interference due to tunnel walls and the model support 
sting at subsonic speeds in slotted-wall wind tunnels. The program is derived from, and uses many elements 
of STIPPAN, a slotted tunnel interference prediction program described in refs. 1 and 2. These programs 
utilize a high order panel method technology as developed by Thomas (ref. 3) augmented by special features 
for representation of slotted-wall tunnel test sections. In PANCOR, use is made of pressures measured on 
the slotted walls between slots to control the flow through the slots and to determine the velocity entering 
the test section. One or more additional wall pressures measured just downstream of the slotted region 
control the velocity entering the solid-wall diffuser. Applications of PANCOR are illustrated in ref. 4. 

The input data are organized in four groups, three of which can be saved for use in successive cases 
and only those data elements most likely to change from one test point to the next in a wind tunnel 
test program need be read for each case. Printed program output includes interference corrections at the 
model, consisting of a longitudinal buoyancy correction to drag coefficient and interference increments of 
Mach number and flow angularity averaged over the wing and over the horizontal tail if one exists. These 
corrections are evaluated both including and excluding sting interference. Standard printed output also 
describes the flow properties at the control points used for problem solution and at a set of field points 
arbitrarily defined by the user. The output properties include the components of tunnel interference 
perturbation velocity. The interference perturbation is defined as the total perturbation less that induced 
by the singularities representing the test model. Output data files accessible to computer plotting utilities 
are also provided. 

The present paper presents a description of the types of panel singularities, networks, and boundary 
conditions provided, gives guidelines for assembling the interference assessment problem from these 
elements, and describes the input data requirements and output format in detail. 

Symbols 

D constant in right hand side of boundary condition 

c wing chord at local spanwise station 

Cp pressure coefficient 

d wall slot spaicing 

rn,n index of m-line or n-line in panel network 

n unit normal vector of panel 

S source line strength 

t local wing thickness 

*TE wing thickness at local section trailing edge 
u longitudinal perturbation velocity 

^oo reference velocity vector 

Vp perturbation velocity vector 



x, y, z cartesian coordinates 

Xi x coordinate relative to mid-chord point of local wing section 
7 local lifting vorticity on wing or tail 

/x doublet strength on panel 

a source sheet strength on panel 

$ perturbation potential in domain exterior to computational tunnel domain 

Network, Panel and Boundary Condition Options 

Network Properties 

Many of the singularity panel networks and boundary conditions used in PANCOR are identical with 
those described in ref. 2 for the STIPPAN program. For completeness, and because some differences exist 
in source networks and boundary conditions, the elements used in PANCOR are described herein in detail. 

In the basic panel method, the bounding surface of a solution domain is divided into one or more 
networks, each of which may be subdivided into panels. The geometry of each network is identified with 
the intersections of a set of m-lines with each of a set of n-lines. The input integers ML and NL define 
the number of lines in each set and the index m or n identifies a specific member of the set of m-lines or 
n-lines respectively. Panel corners are located at the line intersections and the boundaries of each single 
panel are formed by the straight line segments connecting the intersections of two adjacent m-lines with 
two adjacent n-lines. While the panels are quadrilateral in general, adjacent lines in either set are allowed 
to merge forming a triangular panel. 

Within each network, the ML*NL line intersections (defining points) and the (ML- 1)*(NL- 1) panels 
are indexed with the inner index m varying most rapidly. The panel index is used also as the index for 
other properties uniquely associated with each panel. The corners and edges of each network are also 
identified in order. The in and n values for eacli network corner are: 


corner 

m 

n 

i 

i 
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NL 
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ML 

NL 
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ML 
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and the edges are identified by: 


edge 
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NL 
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ML 
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A unit normal vector is calculated for each panel as the unit vector in the direction of the cross 
product of a panel width vector in the n-advancing direction into a panel width vector in the m-advancing 
direction. The sense of the unit normal vector is sensitive, therefore, to the order in which panel defining 
point coordinates are read into the network index grid. 

Certain properties common to all panels within a network are input as network properties. These 
include the type and order of singularity distributions over each panel and the form of boundary condition 
imposed at the control points. 
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Doublet Networks 

The form of doublet strength distribution over each panel is specified for each network by the value of 
IDT. The options available are given in the following table. 


IDT 

Order of 
distribution 

Number of 
unknowns 

0 

none 


1 

constant 

(ML-l)*(NL-l) 

2 

bilinear 

(ML-1)*(NL-1) 

3 

biquadratic 

(ML+l)*(NL+l) 


Control points (points at which boundary conditions are imposed) are located at the panel center points 
for all nonzero doublet distributions. For the biquadratic distribution, additional control points are placed 
at the network corners and along the network edges at panel boundary midpoints. For IDT > 0 the set 
of unknowns for the network is the set of doublet strengths at coutrol point locations. For IDT > 1, 
each coefficient of a higher order term in the bilinear or biquadratic equation for the panel singularity 
distribution is established as a linear combination of the singularity strengths at neighboring control 
points, the combining coefficients having been determined from a weighted minimization process. 

Source Networks 

The options available for panel source strength distributions are given in the following table and are 
selected by the value of 1ST. 


1ST 

Order of 
distribution 

Number of 
unknowns 

Restrictions 

0 

none 


lt>T > o 

1 

constant 

(ML-l)*(NL-l) a 


2 

bilinear 

(ML-l)*(NL-l) a 


6 

bilinear 

NL-1 

ML— 1 a factors specified 

9 

bilinear 
+ lines 

(ML-2)*(NL-2) 

IDEFN=0, IDT=0 

-1 

constant 

1 (free) 

ML=NL=2, IDT=0 or 1 

-2 

bilinear 

3 (free) 

ML=NL=2, IDT=0 or 1 


a for 1DT=0 


Values of 1 and 2 provide general purpose source panel networks with the source strength quantified at 
the panel centers. For 1ST = 2, the bilinear coefficients are established by the same procedure used for 
the bilinear doublet distribution. The set of panel center source strengths is the set of unknowns for the 
network only if IDT = 0 is specified for the same network. If both IDT and 1ST are greater then zero, the 
set of unknowns is determined by IDT and the panel center source strengths are calculated from: 

<r = -Voo • n + B ( 1 ) 

where B is a constant entered into the right hand side of the matrix equation and is zero by default unless 
an option to read a set of B values for the network is selected. 

The 1ST = 6 network is a special bilinear source panel network developed for use in the transition 
region at the downstream end of wall slots where some slotled-wall tunnels have reentry flaps. Control 
points are located only at the center of the last panel in each m row; i.e. the panels between m-lines ML— 1 
and ML. The set of unknowns is the set of source strengths at the centers of these last panels. The source 
strengths of all other panels are established by use of a set of strength facotrs specified by the user. There 
are ML-1 such factors and the ratio of the strength of any panel to that of the last panel in the same 
m-row will be the ratio of the value of the corresponding strength factor to that of the last factor. 
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A discrete slot representation of the slotted wall is provided by the 1ST— 9 network. The discrete slots 
are represented by piecewise linear line source distributions along each n-line except those at network edges 
2 and 4 and are quantified at all interior panel corners. The set of unknowns consists of the line source 
strengths at all interior panel corners. The line source strengths at network edges 1 and 3 are zero by 
definition. Control points are also located at the Ulterior panel corners. Bilinear panel source distributions 
are also used but their strengths and gradients are determined by the line source strengths. IDT must be 
specified as zero but the network must be superimposed on another network having IDT> 0 and having 
panel boundaries in the slotted wall region winch coincide with those in the IST= 9 network. 

1ST values of -1 and -2 may be used only with a single panel network (ML-NL— 2) and offer the means 
of inplemcnting special problem closure conditions. The unknowns consist of the one or three coefficients 
of the constant or bilinear source distribution over the panel. The location and other properties of the one 
or three control points are free to be established directly by input specifications. If IDT— I is specified 
along with a negative 1ST, the doublet strength unknown and the panel center control point are appended 
to the unknowns and control points established by the negative 1ST. 

Boundary Conditions 

The types of boundary condition provided are listed in the following table. Although the boundary 
condition type imposed at each control point is ultimately stored in the array named IBCT, the input 
specification is in a different form. For each network, the type used for all panel center control points 
is given as IBCIP, that used at each edge of the network is taken in order from the four input values of 
IBCEP, and at each corner from the four input values of IBCCP. The boundary condition type for each 
free boundary point provided by negative 1ST, and for each redefined control point is input directly into 
IBCT. 


IBCT 

Description 

Restrictions 

0 

none 

IBCIP, IDT 0, 1IUISI- 1 

1 

Cp specified 


3 

ii 

f 

ta 


4 

0 


5 

/x — U (<r = G if JDT= U) 


6 

slot (lux smoothness 

IBCEP, IDT=0, see text 

12 

dn/Uy = 0 

IBCEP or IBCCP, 1DT> 0 


Boundary condition type 0, if specified as IBCIP, results in the elimination of all control points for the 
network and provides for the a priori specification of source panel strengths in a network with constant or 
bilinear source panels. In this case, IDT must be specified as zero. 

Boundary condition type 1 provides for specification of pressure coefficient at the control point locations 
on the positive normal side of the panels. When this boundary condition is used, the solution is updated 
iteratively to satisfy the exact nonlinear expression for pressure coefficient. Values of prescribed pressure 
coefficient may be read into the right hand side array D. Alternatively, a feature allowing the control 
point properties to be redefined may be used. With this feature, the location, boundary condition type, 
and right hand side constant of any specified control point(s) can be altered from those created from the 
network properties. This feature is particularly useful for pressure coefficient specification to allow control 
point locations to coincide with pressure orifice locations. 

Boundary condition type 3 provides for the direct specification of the normal component of total velocity 
(Neumann boundary condition) on the positive normal side of the panel. The constant B (default=0) is 
the specified value of normal velocity. If 1ST and IDT are both greater than zero, the same value of B is 
used in calculating the panel source strength by eqn. (1). The normal component of perturbation velocity 
on the opposite side of the panel is thereby set to zero regardless of the value of B. 

Boundary condition type 4 provides an indirect means of imposing the Neumann condition. The 
condition imposed directly is that the perturbation potential on the negative normal side of the panel 
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be zero. Again, if 1ST and IDT ?xre both greater than zero, the source strength calculated by eqn. (l) 
produces a normal component of total velocity on the positive normal panel surface approximately equal 
to B . 

Boundary condition types 5 and 12 are special conditions imposed directly on local doublet distribution 
properties rather than on flow properties. Boundary condition type 5 provides a convenient means of 
constraining the otherwise free constant of integration in the value of potential if all other boundary 
conditions are imposed on velocity rather than potential. Boundary condition type 12 can be imposed at 
a doublet network edge lying in a plane of symmetry to improve the solution continuity across the plane 
of symmetry. 

Boundary condition type 6 must be specified as IBCIP. It is provided to suppress an abrupt onset of 
line source strength S at the beginning of a slot line by requiring that the line source gradient dS/dx is 
equal on the first two line segments. On an individual source line, it is formulated as 


x[-x 0 x 2 -x 0 

where the subscripts 0, 1, and 2 denote conditions at the slot intersection with m-lines 1, 2, and 3 
respectively in the 1ST = 9 networks. This boundary condition is not intended for specification on a 
network with 1ST = 9, but may be specified in place of local constraints at control points on networks 
having 1ST — 1 or 6. If the total number of type 6 constraints is equal to the total number of slots in 
all 1ST = 9 networks, eqn. (2) is imposed independently on each slot. Alternatively, if the total number 
of type 6 constraints is equal to the number of 1ST = 9 networks, eqn. (2) is summed over all slots in 
each network. If only one type 6 constraint is specified, eqn. (2) is summed over all slots in all 1ST = 9 
networks. The logic underlying the use of this constraint is discussed in the following section. 

All control points, or constraints, are indexed within each network as shown in the following table. 


Singularity type 

Local index 

n-range 

m-range 

JDT or 131'= 1 or 2 


n=l,(NL-l) 

in— i ,(M L — 1) 

IDT=3 

(n-l)*(ML+l)+m 

n=l,(NL+l) 

m=l,(ML+l) 

1ST=6 

n 

u=l, (NL-1) 


IST=9 

(n-J)*(ML-2)+m 

n=l,(NL— 2) 

m=l,(ML— 2) 


The global control point index is formed by adding the index within the local network to the total number of 
all preceding control points. For this purpose, all control points in networks with boundary condition type 
6 are accumulated first, then those in all other networks are accumulated in the order of their definition 
in the input data file. The order thus established for the global control point index is used later for the 
order of rows in the matrix equation. 

Smoothing 

PANCOR provides a capability for altering the primary coefficient matrix to introduce solution 
smoothing within any multi-panel network. The user may specify smoothing factors in the in- and n- 
line directions of any such network which causes the solution distribution to be smoothed in the specified 
direction while allowing some error in satisfying the corresponding boundary conditions. This feature is 
useful to allow solution of an otherwise divergent problem or to improve the regularity of the solution of 
an ill-conditioned problem. 

Structuring the Assessment Problem 

Program PANCOR provides a number of features, or building blocks, which can be linked together to 
define and then solve a wind-tunnel interference assessment problem. Some of these building blocks were 
described in the previous section; in this section, their appropriate assembly will be discussed. Although 
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numerous variants are possible, the two variations described herein are those which have received the most 
emphasis in development of the PANCOR program. 

The solution domain, occupying a rectangular parallelepiped, is set up as a fully bounded domain in 
which the potential flow field represents the equivalent inviscid flow in a portion of a wind tunnel including 
the test section, a constant area solid walled entrance duct, and the initial part of the solid walled diffuser 
just downstream of the test section. Outside of this domain, an unperturbed flow is presumed, having a 
velocity of unit magnitude in the direction of the tunnel axis. This unit velocity is taken as the reference 
velocity from which perturbations are expressed in both the solution and the outer domains. The solution 
domain is bounded by networks of doublet panels which permit a potential jump across the boundary. Type 
4 boundary conditions are specified with these panels to fix the outer flow in an essentially unperturbed 
state. The potential in the solution domain is made continuous with that in the outer How by specifying 
zero strength of a doublet panel on the upstream face of the solution domain. If this face is sufficiently far 
upstream of the test section or other sources of disturbance, it may be represented by a one-panel network 
(NL=ML=2) in which zero doublet strength is achieved by omitting the doublet panel (IDT=0). 

Boundary conditions are imposed on the interior flow by superimposing one or more source networks 
on the boundary doublet networks. If Cp were known over the entire bounding surface, source networks 
with the type I boundary condition could be used to reproduce the tunnel flow. The LIN COR program 
(ref. 5) presumes that this kind of boundary information is known. Program PANCOR, however, allows 
mixed boundary conditions, thereby permitting the use of the Neumann condition where the equivalent 
inviscid normal velocity component at the boundary is known with sufficient accuracy, and requiring the 
use of measured pressures only as needed to establish flow rates normal to the boundary in those regions 
where the flow processes are not known with confidence. Neumann conditions are imposed on solid wall 
regions of the tunnel by specifying boundary condition type 0 for the source panel network and specifying 
the right hand side constant D at each panel center as the equivalent wall slope. The Neumann condition 
in regions with zero wall slope can be imposed simply by omitting any source panels in those regions. 

Recourse is made to measured pressures to define the flux distribution through the wall slots and 
through the reentry region at the downstream end of the slots, and the longitudinal velocity perturbation 
of the flow entering the upstream end of the solution domain. To control the flux distribution through the 
slots, pressures measured on the flat surfaces between slots are specified in coefficient form. The number 
of longitudinal rows of such pressures must equal the number of slots. The slotted wall is represented in 
PANCOR by superimposing the discrete slot source network, IST=9, on the doublet panels in the slotted 
wall regions. All of the control points in the IST=9 network are then redefined to be located laterally in a 
measured pressure row, and longitudinally, some fraction of panel length downstream of the corresponding 
line source quantifying point. The magnitude of the longitudinal control point shift is not critical although 
solution stability is best if the control points are about one-half panel length downstream of the line source 
quantifying points. 

At the downstream end of the slots, some slolted-wall tunnels have a region of complex geometry 
involving reentry flaps and possibly a discontinuity in wall contour. It is difficult to state just how such 
a region should be modelled or just how much use should be made of measured pressures in this region. 
With PANCOR, it is possible to model such a region with a smooth transition from discrete slot flux to 
smoothly distributed wall flux. The IST=9 network should extend to the downstream end of the slots but 
if the next to last in-line is located at the leading edge of the reentry flap (or transition region) the line 
source strength will decrease linearly from that at the last quantifying point to zero at the slot end. An 
1ST=6 network is then superimposed on this region with the shape of the longitudinal source distribution 
specified to increase smoothly (perhaps linearly) to that of the downstream panel. Wall pressures measured 
in the immediate vicinity of the slot ends may be specified with a type 1 boundary condition to control 
the strengths of the longitudinal rows of panels in the IST=6 network. It is preferable to divide the IST=6 
network laterally so that each slot terminates in a different longitudinal row of panels. Care should be 
exercised, however, to avoid placing a panel boundary very close to a pressure control point. 

Figure 1 illustrates the superposition of the several kinds of panel networks used to represent a slotted 
wall. As illustrated in this figure, large source panels with specified strength may be added to represent 
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regions of constant wall slope. The combination of the IST=9 discrete slot network and a IDT=3 doublet 
network produces essentially zero normal velocity everywhere on the slotted wall except at the slots 
themselves, regardless of the slot flux. Thus, the effective slope of a slotted wall may be controlled by 
adding specified strength source panels just as on a solid wall. 

One more property of slotted walls must be recognized in the PANCOR problem formulation, that is, 
the flow must make a smooth transition from the solid wall nozzle to the slotted wall region without an 
abrupt onset of flux through the slots. This property, which is not unlike the Kutta condition at an airfoil 
trailing edge, is demonstrated by the characteristic mode of tunnel/plenum interaction derived in ref. 1. 
This property may be satisfied in PANCOR by using a type G boundary condition to control the source 
strength of a panel on the upstream face of the solution domain, thereby allowing a longitudinal velocity 
perturbation to exist in the flow approaching the slot origin to compensate for a difference in definition of 
the reference pressure used in forming pressure coefficients between the experimental data reduction and 
the PANCOR computation. On the downstream face, a one-panel network with 1DT= 1 and IST= — 1 
may be used to impose the unperturbed outer flow constraint, bounday condition type 4, on the centers 
of both the downstream and upstream faces of the computational domain. 

In the formulation just described, the slot flux smoothness constraint is imposed only in an average 
sense over all slots. Minor inaccuracies in the measured wall pressures can cause erratic flux distributions 
on individual slots. Also, the pressures used to control the 1ST=6 source strengths must be near the slot 
terminations at the entrance to the solid wall diffuser. Because of the complex flow phenomena occurring 
here, it is anticipated that the spacial variation of pressure in this region of the actual tunnel might be 
large and might be poorly reproduced in the computational flow. Thus, adaptation of this PANCOR 
formulation to a given tunnel might require careful tailoring of the locations of both the wind tunnel 
orifices and the PANCOR control points. 

An alternative formulation which reduces some of these problems can be used. In this alternative, slot 
flux smoothness constraints are specified for each control point in the IST=6 networks at the downstream 
end of the slotted walls. Measured pressures are not used at these points. Instead, a single wall pressure, 
measured at a location on the solid wall diffuser where the actual and computational pressures might be 
matched more reliably, is used to control the source panel strength on cither the upstream or downstream 
face of the computational domain. This alternative formulation is illustrated in fig. 1 and is used in the 
sample case described in Appendix D. It probably is preferable in most cases to that using measured 
pressures as boundary conditions for the 1ST=6 networks. 

The Input Data File 

The input data file is a sequence of card image records. The following subsection entitled Input File 
Records describes in order all of the record types which are or might be required. The input variables, 
format, and repetition requirements for each record type are given and the conditional requirements are 
noted in comments. If the number of records is indicated as a fraction, the fraction should be rounded 
upward to the next integer. 

A following subsection lists the definitions of all input variables in the order encountered in the input 
file description. Note that record type 1 is read only as the first line of the input file; the value of CASID, 
if greater than zero, is automatically incremented by unity upon completion of each case. Input for the 
second and succeeding cases begins with the MOVE array in record type 2. After the first two record 
types, the input data are divided into four groups, the first three of which may or may not be included 
depending on the entries in the MOVE array. All groups must be included for the first case of a file, and 
the MOVE array and the case-dependent group must be included for each case. 

The appendices A, B, and C give more detailed descriptions of the MOVE array, the network linkage 
provision, and the test model representation, respectively. A sample case is described in Appendix D 
which shows a complete input data file. 
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Input File Records 


Record 

No. or 



Type 

Records 

Format 

Variables or comments 

I 

1 

F10.0 

CAS1D - First line only, do not repeat for multiple case input files. 

2 

1 

414 

MOVE(l-4) 




The panel definition group, record types 3 through 15, is to be 
included only if MOVE(4)— 0. 

3 

1 

2I5.2F10.0 

1SYM, NNET, XROT, ZROT 

4 

NNET 

1714 

NL, ML, NLIl, IDEFN, 1ST, IDT, IRHS1, NCPR, IBC1P, 
lBCEP(l-4), lDCCP(l-4) 

5 

1 

I10,2F 10.0 

ITMX, CNVU, RFU 

6 

NNET/8 

8F10.0 

FSN(l) 

7 

NNET/8 

8 F 10.0 

FSM(I) 




For each network, taken in order from 1 to NNET, include either 
a group of type 8 through 10 records if JDEFN=0, or the required 
number of type 11 records if 1DEFN = 1. 

8 

l 

3I5,F 10.0 

NDIR, MDJR, NORM, XNO 

9 

NL/8 

8F10.0 

XND(I) 

10 

ML/8 

8F10.0 

XMD(I) 

11 

NL*ML 

3F10.0 

PDEFP(I) 

Include a record type 12 only for each network having IST=6 and 
IRI1S1— 0. 

12 

(ML-l)/8 

8F10.0 

SFM(I) 




Include the required number of type 13 records for each network 
having 1RHS1=1. 

13 

(NL-1)* 

8F10.0 

BCENP(l) 


(ML-l)/8 


NFBP is total number of free boundary points. Count 1 for each 
1ST= —1 network and 3 for each IST= —2 network. 

14 

NFBP 

215.6F10.0 

IPOINT, 1BCT, PCONP(l-3), SCONP(l-3) 

15 

NCPR 

21 5,3 FI 0.0 

1CONP, IBCT, PCONP(I-3) 




The field survey group, record types 16 through 18, is to be 
included only if MOVE(2)=l. 

16 

1 

315 

NROW, IDAT, NDAT 




Include the required number of type 17 records only if IDAT>0. 

17 

N DAT/8 

8F10.0 

DAT(I) 

18 

NROW 

215, 6F 10.0 

IPCI1, NPROW, XI, Yl, Zl, X2, Y2, Z2 


The model and sting group, record types 19 through 34, is to 
be included only if MOVE(3)-l. Note that PANCOR does not 
provide for model roll attitudes other than the conventional 
upright attitude. An inverted model (180° roll) may be simulated 
by reversing the sign of the variables underlined below in record 
types 20, 21, 23, 27 and 37. If the sting is also inverted, the values 
of ZSS in record type 32 should also have reversed sign. 
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19 

1 

515 

NBS, NWS, NTS, NSEB, ISEB 

20 

1 

5F10.0 

XMROT. ZMROT. DTI1ET. ZWING. ZTAIL 

21 

1 

5F10.0 

XMREF. ZMREF. SREF. CHEF, OA2MI 




Include types 22 through 25 only if NBS> 2, omitting type 25 if 
NSEB— 0. 

22 

NBS/8 

8F10.0 

XBS(I) 

23 

NBS/8 

8F10.0 

ZBS(U 

24 

(NBS-1J/8 

8FI0.0 

QBV(I) 

25 

NSEB/8 

8F10.0 

WKW(I) 

Include types 26 and 27 only if NWS>2. 

26 

NWS 

7F10.0 

YWG, XCW, GW, QS0, QS1, QS2, QS3 

27 

NWS 

8F10.0 

QG01. ASH. QG21, QG31, QG02. QG 12. QG22. QG32 
Include types 28 and 29 only if NTS> 2. 

28 

NTS 

7F10.0 

YTS, XCT, CT, QS0T, QS1T, QS2T, QS3T 

29 

NTS 

4F10.0 

QG0T, QG1T, QG2T, QG3T 




Type 30 required if MOVE(3)- 1. 

30 

1 

315 

NST, NSEP, ISEP 




Include types 31 through 34 only if NST> 2, omitting type 34 if 
NSEP=0. 

31 

NST/8 

8F10.0 

XSS(l) 

32 

NST/8 

8F10.0 

ZSS(I) 

33 

(NST-l)/8 

8F10.0 

QSV(I) 

34 

NSEP/8 

8F10.0 

WKS(I) 




The case-dependent group, record types 35 through 40 is to be 
included for all cases, noting the exceptions given. 

35 

1 

10A8 

TITLE 

36 

1 

4 F 10.0 

TEST, TRUN, TPNT, TMACH 

37 

l 

5 F 10.0 

AMIIEF. TI1ETS. CLIFT. CDRAG. CMOM 


38 


39 


(NL-1)* 

(ML-l)/8 


8F10.0 


At this point, include a record of type 12 for each network having 
IST=6 and 11U1SI=3. 

Include the required number of type 38 records for each network 
having IRIISI=2. 

BCENP(l) 


40 NCPRT/8 


Include a record type 39 for each network liaving NFBP> 0. 
8FI0.0 BCONP 

NCPRT is summation of NCPR over all networks. 

8F10.0 BCONP(I) 
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Definition of Input Variables 


CAS1D 

MOVE 

1SYM 


Case identification number to be recorded on SIF file. Should have 
integer value for proper SIF file usage. Value of 0.0 causes no SIF file 
to be written. 

A 4-eleinent integer array for program control (see Appendix A). 

Symmetry flag. 

= 0 for no symmetry. 

= 1 for symmetry about y=0 plane. Model body and sting are not reflected. 


NNET 

XROT, ZROT 
NL, ML 
NLR 
IDEFN 


1ST 

IDT 


Total number of networks. 

x- and z-coordinates of center of rotation of sting support system. 
Number of n-lines and m-lines for a network. 

Receiving network number for linked output (see Appendix 13). 

Panel geometry input flag. 

= 0 for simplified orthogonal network input form. 

= 1 for input listing of panel corner point coordinates. 

Source distribution type. 

Doublet distribution type. 


1RIISI 

= 0 
= 1 

= 2 
= 3 

NCPR 
I POIP 
IBCEP(l) 

IBCCP(I) 

1TMX 

CNVU 

RFU 

FSN 

FSM 

NDIR, MDIR, NORM 

XNO 

XND(I) 

XMD(I) 

PDEFP(l) 

SFM(I) 

BCENP(I) 


Right hand side constant input flag, 
for B = 0. 

for B read from record type 13 in paneling difinition group, 
for B read from record type 38 in case-dependent data group, 
for SFM read from record type 12 in case-dependent data group. 

Number of redefined control points in network. 

Boundary condition type at network interior control points. 

Boundary condition type at control points on network edge I. 

Boundary condition type at control points at network corner I. 

Maximum number of iterations allowed for satisfaction of exact pressure 
coefficient boundary conditions. 

Convergence criterion for maximum residual of nonlinear terms in u as 
function of Cp. 

Relaxation factor for update of exact Cp boundary conditions. 
Smoothing factor in direction of varying n. 

Smoothing factor in direction of varying ni. 

Coordinate direction of varying n-index, varying m-index, and network 
normal respectively. Use 1, 2, 3 for x, y, z. 

Value of NORM coordinate at network plane. 

Array of n-line coordinates in NDIR direction (1=1 to NL). 

Array of in-line coordinates in MDIR direction (1=1 to ML). 

Coordinates of panel corner points in network (1 = 1, 2, 3 for x, y, z). 

Array of panel center source strength factors for all longitudinal panel 
rows in 1ST=6 network. 

Array of right hand side constant B for all panel center points in 
network. 
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1POINT 


Panel index number to define panel normal recession direction of free 
boundary point. 

IBCT 


Boundary condition type. 

PCONP(I) 


Control point coordinates (1=1, 2, 3 for x, y, z). 

SCONP(I) 


Components of unit normal vector used for boundary condition type 3 
(1=1, 2, 3 for x, y, z). 

ICONP 


Local control point index within network. 

NROW 


Number of straight rows of (low survey points. 

IDA'I' 

= 0 

Uniformity (lag for flow survey point spacing, 
for uniform point spacing in all rows. 


> 0 

use 1, 2, or 3 for x, y, or z coordinate values, respectively, given in DAT. 

NDAT 


Number of non-uniformly spaced coordinate values given in DAT. 

DAT(l) 


Array of coordinate values of non-uniformly spaced points projected on 
the axis indicated by ID AT. Include only those points lying between the 
first and last points of the row. 

IPCII 

= 0 

Coordinate system and uniformity flag for flow survey points, 
for uniformly spaced points defined in tunnel coordinates. 


= 1 

for uniformly spaced points defined in model coordinates. 


= 2 

for non-uniform points defined in tunnel coordinates. 


= 3 

for non-uniform points defined in model coordinates. 

NPROW 


Number of points in a row of flow survey points. If IPCII=2 or 3, 
NPROW must equal NDAT+2. 

XI, Yl, Z1 


Coordinates of first point in a row of flow survey points. 

X2, Y2, Z2 


Coordinates of last point in a row of flow survey points. 

NBS 


Number of body stations. 

NWS 


Number of wing stations. 

NTS 


Number of tail stations. 

NSEB 


Number of consecutive body segments having separated flow. 

ISEB 


Index of initial body segment having separated flow. 

XMROT, ZMROT 


x- and z-coordiuates of center of rotation in model coordinate system. 

DTIIET 


Pitch angle of model coordinate system relative to sting, degrees. 

ZWING, ZTAIL 


z-coordinate of wing reference plane or tail reference plane respectively 
in model coordinate system. 

XMREF, ZMREF 


x- and z-coordinates of moment reference point in model coordinate 
system. 

SREF 


Model reference area. 

GREF 


Model reference chord. 

DA2M1 


Angle-of-attack change from that corresponding to first lift coefficient to 
that corresponding to second lift coefficient, degrees. 

XBS(I) 


Array of x-coordinates of body stations in order of increasing x in 
model coordinate system (1=1 to NBS). 

ZBS(I) 


Array of z-coordinates of body stations in model coordinate system 
(1=1 to NBS). 

QBV(I) 


Volume of body segment between stations I and I I I (1=1 to NBS-l). 



WKW(I) 

YWG(I) 

XCW(I) 

CW(I) 

QSO(I), QSI(I) 
QS2(I), QS3(I) 

QG01(I), QG1 1(1) 
QG21(I), QG31(1) 

QG02(I), QG 12(1) 
QG22(I), QG32(I) 

YTL(I) 

XCT(l) 

CT(I) 

QSOT(I), QS1T(1) 
QS2T(1), QS3T(I) 

QGOT(I), QG1T(I) 
QG2T(1), QG3T(I) 

NST 

NSEP 

ISEP 

XSS(I) 

ZSS(l) 

QSV(i) 

VVKS(I) 

TITLE 

TEST, TRUN, TP NT 

TMAC1I 

AMREF 

TIIETS 

CLIFT 

CD RAG 

CMOM 

RCONP(I) 


Array of wake widths behind body segments with separated How (1=1 
to NS ED). 

Array of y-coordinates of wing stations in order of increasing y (1=1 to 
NWS). 

Array of x-coordinates of wing station inid-cljord points in model 
coordinate system (1=1 to NWS). 

Array of local wing chords (1=1 to NWS). 

Coefficients of multipole representation of wing section thickness 
distribution at wing station J (see Appendix C). 

Coefficients of inultipole representation of wing chordwise circulation 
distribution at first lift coefficient at wing station 1 (see Appendix C). 

Coefficients of inultipole representation of wing chordwise circulation 
distribution at second lift coefficient at wing station 1 (see Appendix C). 

Array of tail station y-coordinates in increasing order (1=1 to NTS). 

Array of x-coordinates of tail station mid-chord points in model coordi- 
nate system (1=1 to NTS). 

Array of local tail chords (1=1 to NTS). 

Coefficients of mullipole representation of tail section thickness 
distribution at tail station I (see Appendix C). 

Coefficients of mullipole representation of tail chordwise circulation 
distribution at tail station 1 (see Appendix C). 

Number of sting stations. 

Number of consecutive sting segments having separated flow. 

Index of initial sting segment having separated flow. 

Array of sting station x-coordinates at TIIETS=0 relative to center of 
rotation in increasing order (1=1 to NST). 

Array of sting station z-coordinates at THETS=0 relative to center of 
rotation (1=1 to NTS). 

Volume of sting segment between stations I and 1 (-1 (1^=1 to NST- J). 

Separated How wake width behind sting segment between stations 
I I ISEP -1 and H ISEP (1=1 to NSEP). 

80 character case identification label. 

Ileirarchic case identifiers to be passed to TAPE3. 

Tunnel Mach number passed to TAPE3. 

Reference Mach number for PANCOR solution. 

Pitch angle setting of sting support system, degrees. 

Model lift coefficient. 

Model drag coefficient. 

Model pitching moment coefficient. 

Array of right hand side constant B. 
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Array Size Limitations 

Program dimensions limit the maximum size of certain input variable arrays and combinations thereof 
as follows: 


Quantity 

Name 

Maximum 

Number of networks 

NNET 

20 

Number of panel defining points 


1200 

Total number of unknowns 


600 

Number of redefined control points 


300 

NL for each IST—9 network 


8 

NL for each IDEFN=0 network 


20 

ML for each 1DEFN=0 network 


25 

Number of body stations 

NBS 

50 

Number of wing stations 

NWS 

20 

Number of tail stations 

NTS 

20 

Number of sting stations 

NST 

70 

Number of field survey points 


600 


In addition, the number of stored influence coefficients must not exceed 1,600,000. The actual number 
written to storage is given by the output variable MLWD (see fig. 3b concluded). 

Program Organization and Computer Interface 

Program PANCOR is written in FORTRAN 77 with minor exceptions and is suitable for automatic 
vectorization with the CDG' FTN200 Cycle 670 compiler and the CDC VSOS Version 2.3 operating system. 

The program is made up of 24 code modules linked by the calling paths shown in fig. 2. The MAIN 
program is a simple executive routine which calls the seven major groups of subroutines in sequence 
and records the CPU time utilized in each. The INPUT group of subroutines reads all input data and 
deals with panel and control point geometry. The SFIT group produces data relating the higher order 
coefficients of the singularity distribution on each panel to the singularity strengths at neighboring control 
points. These data are written panel by panel to a scratch file identified sis TAPE1 which is accessed 
in the MATA, PCOUT, and FIELD subroutines. The RECESS routine performs final adjustments to 
the control point locations. The next group of routines generates much of the problem forcing data 
accumulated on the right hand side of the matrix equation including those representing the test model and 
support sting. In the MATA group, the aerodynamic influence coefficients for all panels and source lines 
are calculated, accumulated according to the higher order singularity coefficients, stored for subsequent 
use and assembled according to boundary condition type into the primary coefficient matrix. All of these 
operations are performed in a single network loop to minimize data paging into and out of core memory. 
The basic influence coefficient storage array has a dimension of 1,600,000 and is accessed again in subroutine 
PCOUT. A smoothing matrix is created if called for and is summed with the primary coefficient matrix. 

The matrix equation is solved in the MATSOL group. The solution is iterated to update the exact 
pressure coefficient boundary conditions. The gaussian elimination subroutine used is a Langley math 
library routine which factors the linear coefficient matrix only for the first solution and simply performs 
a back substitution for subsequent iterations. The PCOUT subroutine prepares the basic solution output 
data at the panel network control points. The FIELD group uses the singularity strengths defined by 
the solution to produce a How survey at new field points, either calculated by the program for evaluating 
model corrections, or arbitrarily specified by the user. This requires the calculation of new aerodynamic 
influence coefficients but data handling is minimized by accumulating the results directly in the output 
arrays. 
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Five files are opened by the program. TAP El is a binary file which, as was previously noted, is a 
scratch file written and read by the program. The remaining files are coded files. The input data file is 
identified as TAPES, and TAPE6 is the file used to format the printed output. TAPE2 and TAPE3 are 
special output files used to convey solution data to plotting utilities. They are written in the format of a 
Transferable Output ASCII Data (TOAD) file as described in ref. G. 


Toad File Output 

TAPE2 contains solution and field survey results in the TOAD file format and is intended for 
postprocessing by one of the graphics utility programs in the Langley Research Center system of data 
processing utilities. The names in the TAPE2 file LA11EL record are: 

CASE, NETWORK, ROW, POINT Ileirarchic data identifiers. 

XLOC, YLOC, ZLOC Point coordinates. 

VXTOT, VYTOT, VZTOT Components of total velocity. 

VXINT, VYINT, VZINT Components of interference velocity. 

PCOEF Value of pressure coefficient. 

PlIJ Perturbation potential. 

SDX Longitudinal integral of panel source strength. 

DELM, DALPII Interference increments in Mach number and angle of 

attack . 

DELMN, DALPHN Interference increments omitting sting interference. 

The value in CASE starts with CASID from record 1 of the input data file and is incremented by unity 
for subsequent cases in a multi-case run. NETWORK values from 1 to NNET give data at panel center 
control points at ROW values from 1 to NL-1 and POINT vulues from 1 to ML-1. Generating networks 
in a linkage set are omitted unless 1ST=9. Data given for an 1ST=9 network are at slot control points 
at ROW values from 1 to NL-2 and interference velocity components are not given; instead, VXINT 
contains the x location of the line segment end where S is quantified and VZINT contains the line source 
strength normalized to an equivalent homogeneous normal velocity, —S/2d. Data from field survey rows 
are identified by a NETWORK value of NNET+1 and are given for ROW values from 1 to NROW and 
POINT values from 1 to NPROW. The TAPE2 output just described is the standard output included if 
MOVE(l) is input as 1 or more. If MOVE(l) is 0 or 1, the TAPE2 output is limited to that from the field 
survey rows. 

TAPE3, also written in TOAD file format, contains the interference corrections at the model together 
with identifying and correlating data. The names given to the data are: 


TEST, RUN, POINT 
MACH, CL 
DELMW, DALPW 
CDB, CDU 
DCMUP 


DCLSC, DC'MSC 

DELMT, DALPT 
DELMWN, DALPWN 
CDDN, CDUN 
DCMUPN 

DCLSCN, DCMSCN 
DELMTN, DALPTN 


Case identifiers. 

Correlating data for plot abscissa. 

Mach and alpha corrections at wing. 

Drag coefficient corrections for buoyancy and upwash. 
Pitching-moment coefficient correction for span wise variation of 
upwash. 

Lift and pitching- moment coefficient corrections for streamline 
curvature. 

Mach and incidence corrections at tail. 

Corrections as above but excluding sting interference. 
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The TAPE3 file is small, containing only a single set of data values for each case in the PANCOR run. 
Because the TOAD file format is editable, however, the TAPE3 files from many PANCOR runs may be 
merged for convenience in making comparative plots. The TAPE2 file, on the other hand, contains many 
data records for each case. If many of these files are to be saved for future analysis, it might be desirable 
to convert tliam to SIF files (see ref. 7) which are binary files requiring only two-thirds the storage space 
of the corresponding TOAD files. 

Printed Output 

The format of the printed output from program PANCOR is illustrated in fig. 3. The output shown 
resulted from the sample case described in Appendix D. Fig. 3, however, includes only enough of the 
output to illustrate the format and identify the terminology used in output headings. Fig. 3a shows the 
standard output which is printed if MOVE(l) is set to one. The additional geometry details printed if 
MOVE(i) is set to 2 or 4 are illustrated in fig. 3b. The more detailed solution results printed if M0VE(1) 
is set to 3 or 4 are illustrated in fig. 3c. The additional solution results include a listing headed “Flow at 
Control Points” in which the results given at points having boundary condition type 4 apply to the exterior 
domain. If MOVE(l) is set to zero, the printed output is reduced to that pertaining to the field surveys, 
headed only by the case number and title and the last line of the iteration history. If a model exists, an 
appropriate set of field points is generated automatically and used to calculate the data corrections written 
to the TAPE3 TOAD file as described in the preceding section. These corrections are printed as part of 
the field point output. If MOVE(l) is set to -1, the printed output consists only of the model corrections 
and the identifying header data. 

A number of input quantities are repeated in the printed output to aid in case identification and input 
verification. These quantities are identified in the output by the same variable names used for the input 
data file and are defined in a preceding section. Additional output quantities are defined below. 


Definition of Output Quantities 


Standard output, MOV E( l)~ I 


NNE 

NDEFP 

NDEFPL 

NCENP 

NCENPL 

NCONP 

NCONPL 

NSSP 

NSSPL 

NDSP 

NDSPL 

1CONP 

QDS 

QCD 

QSS 

QSD 

ITER 

1UMX 


DUMX 


Network index. 

Number oT defining points in network. 

Cumulative number of defining points in preceding networks. 

Number of panels (center points) in network. 

Cumulative number of panels in preceding networks. 

Number of control points in network. 

Last global control point index preceding current network. 

Number of source singularity panels in network. 

Cumulative number of source panels in preceding networks. 

Number of doublet singularity panels in network. 

Cumulative number of doublet panels in preceding networks. 

Global control point index. 

Strength of point source in model body representation. 

Strength of line doublet segment in model body representation. 
Strength of point source in sting representation. 

Strength of line doublet segment in sling representation. 

Iteration step number. 

Control point having largest change in u in current iteration, identi- 
fied by position in ordered consecutive array of all control points with 
type 1 boundary condition, 
u change in current iteration at IUMX point. 
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U1 


STIME 
X, Y,Z 


VX, VY, VZ 

VX1NT, VYINT, VZINT 
CP 

INTSDX 

S 

piii 

I ROW 
POINT 
DEEM 
DALPII 

TIME 

DTIME 

PETIME 


u perturbation at upstream closure panel. For correct value, the 
upstream closure panel must be defined by the first type 4 record in 
the input data file. 

Cumulative computer CPU time used for matrix equation solution. 
Nominal coordinates or control point where results are given, solution 
domain side of panel is implied in panel center output listing, and 
slotted wall output listing. 

Components of total velocity at control point. 

Components of interference velocity at control point. 

Pressure coefficient. 

Integral of panel source strength in x direction from m=l network 
line to present position. 

In slotted wall output listing, line source strength at line source 
quantifying point. 

Perturbation potential at X, Y, Z. 

Row index identifying FIELD survey row. 

Point index in FIELD survey row. 

Increment in Mach number due to tunnel and sting interference. 
Increment in angle of attack (degrees) due to tunnel and sting 
interference. 

CPU time in seconds from start of job. 

CPU time increment used in each group of subroutines. 

CPU time increment used in each group of subroutines expressed as 
percent of CPU time for current case. 


Additional geometry output. M0VE(l)-2 or 4 


IDEFP 
NPB 
X, Y, Z 

XN, YN, ZN 
AREA 
NNNl 
NNN2 

FP(X, Y, Z) 

1SF 
NS 
ND 
1.1 F 


NFP 

MLWD 

LPF 


Global index of defining points. 

Number of panel boundaries. 

Panel center coordinates. (Panel singularity lit output is in incom- 
pressible domain). 

Components of panel unit normal vector. 

Panel area, in incompressible domain. 

Number of neighboring control points used in bilinear singularity fit. 
Number of neighboring control points used in biquadratic singularity 
fit. 

Final coordinates of control points after recessing from panel surfaces 
and network edges. 

Singularity flag indicating source, doublet or both. 

Number of terms in panel source distribution equation. 

Number of terms in panel doublet distribution equation. 

Flag indicating existence of control points exactly coplanar with any 
panel in network. If flag is set, a warning message is issued indicating 
number of control points (including those reflected by symmetry) 
coplanar with that panel. Warning may be ignored if all coplanar 
control points lie outside the pauel boundaries. 

Number of influenced points for influence coefficient computation. 
Cumulative number of words stored in influence coefficient storage 
array. 

Cumulative number of large pages of virtual memory filled by influ- 
ence coefficient array. 
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Additional solution 

NRM 

GENP 

S 

CSX, GSY, GSZ 
I) 

GDX, GDY, GDZ 


output, MOVE(l)=3 or 4 

Number of rows in matrix equation. 
Global center point index. 

Source strength at panel center. 
Components of panel source gradient. 
Doublet strength at panel center. 
Components of panel doublet gradient. 



Appendix A 
The MOVE Array 

MOVE is a 4-element integer array used for program control. The available options for each element 
are listed below. 


MOVE(l) 


Printed output control. 


= -1 

Model correction output only. 


= 0 

Field survey and model correction output only. 


= 1 

Standard output. 


= 2 

Additional printed geometry output. 


= 3 

Additional printed solution output. 


= 4 

Additional printed geometry and solution output. 

MOVfi(2) 


Field survey control. 


= 0 

No field computation. 


= 1 

Read new field specifications. 


= 2 

Use field specifications from previous case. 


= 3 

Same as 2 (reserved for future use). 


= 4 

Used internally to end case after divergent iteration. 


= 5 

Used internally to end job after last case. 

MOVE(3) 


Model and sting control. 


= 0 

No test model. 


= 1 

Read new model and sting specifications. 


= 2 

Use model and sting specifications from previous case. 

MOVE(4) 


Process control. 


= 0 

Read new paneling geometry and execute complete program 


> 1 

Use previous case paneling geometry. 


> 2 

Use previous case singularity fit. 


> 3 

Use previous case aerodynamic influence coefficients. 


= 4 

Use previous case factored matrix. 


Savings in computer resources can be achieved by executing multiple cases in a single job submission 
with judicious use of MOVE(4). Constraints on the allowable use of MOVE(4) are given below. 

MOVE(4)=0 must be used for: 

a. First case in input file. 

b. Any change in panel geometry or other data in input record types 3 through 15. 

MOVE(4)< 1 must be used for any case in which Mach number difTers from the previous case. 

MOVE(4) values of 2 and 3 were used in program development but are of little use with the program in 
its present form. 

MOVE(4)~ 4 requires the least computing time and may be used for all cases not requiring MOVE(4) 
values of 0 or 1. 
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Appendix B 

Network Output Linkage 

The listing of solution output gives the flow characteristics on the positive normal side of the panel at 
the center point of each panel in each network, and optionally includes the panel singularity strengths and 
gradients. If the boundary condition is type 4, the aerodynamic influence coefficients describe the flow 
properties on the opposite side of the panel. The local doublet strength and gradients and source strength 
are then used to transfer the flow potential and velocity components to the positive normal side of the 
panel. Network output linkage provides the capability to combine the source strengths and gradients of 
all coplanar networks within eacli linked group. 

The program assigns a network number NNE in sequence from one to NNET in the same order in 
which the networks are defined by the type 4 records in the input data file. Groundrules for the use of 
linkage are fisted below. 

1. Linkage occurs in groups with one receiver network and one or more generator networks in eacli group. 

2. Doublet panels may exist only in the receiver network and must not exist in the generator networks. 

3. The network number of each generator network in a linkage group must be higher than that of the 
receiver network for that group and lower than that of the receiver network for the next group. Non- 
linked networks may be interspersed at will. 

4. Linkage is invoked by setting NLR for each generator network in one group equal to the receiver network 
number for the same group. NLR for each receiver or noil-linked network should be set to zero. 

f>. AH networks involved in linkage must be defined with 1DEFN-0 and those in each group must have 
identical values of ND1R, MDIlt, NORM and XNOR. (Only flat, orthogonally oriented, coplanar 
networks may be linked together.) 
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Appendix C 

Model and Sting Representation 


Distributed Singularity Model 


The representation of the test model provided in PANCOR is evolved from that used in program 
LINCOR by Rizk and Smithmeyer (ref. 5) and is the same as that used in ST1PPAN (ref. 2), but is extended 
in PANCOR to provide for matching specified values of lift, drag and pitching moment coefficients. The 
basic representation is described more fully in ref. 1. ihe model consists of three components, body, wing 
and tail. Each component i 3 defined by input data given at a specified number of stations. If the number 
of stations given Tor any component is less than two, no further input data is read for that component 
and it makes no contribution to the model perturbation. For convenience of geometry input, separate 


reference coordinate ayylems are used for the inodel and for the stiug. Ihese two systems, together witli 
the wind tunnel coordinate system share a common plane of symmetry at y-0. The center uf sting 
rotation is the origin of the sling coordinate system and is located in the inodel reference system by the 
coordinates XMROT and ZMROT. The angle DT11ET is the pitch orientation of the model reference 
system relative to the sting axis. r l he model and sting are then located in the tunnel by the sting rotation 
center coordinates XROT and ZROT in the tunnel coordinate system and the sting pitch angle THETS. 


The Mode I Body. The body representation makes use of inclined slender body principles in which a point 
source represents a change in cross section area scaled by cosine of angle of attack, and a line doublet 
segment represents the local cross section area scaled by sine of angle of attack. The present program 
applies this concept segment by segment to accommodate an irregular body camber shape. The body 
axis is located in the y=0 plane. Body input data should describe the full body rather than a half body 
because body influence computations are independent of the input value of ISYM. The input quantities 
NBS, XBS, and ZBS give the number and coordinates of stations along the body axis and the volume 
of each segment between stations is input into QBV. A separated wake is presumed to trail from the 
blunt base of the last body segment. The wake displacement is the cross section area uf the last segment 
scaled by cosine of the segment angle of attack. '1. his wake may be eliminated by appending a dummy 
segment having zero volume to the end of the body. Integration of the body with a sting is discussed in 
the subsequent section describing the sting representation. 

Capability is provided to represent additional wake blockage due to flow separation from inclined body 
segments. NSEB is the number of consecutive body segments generating a wake, ISEB is the segment 
number of the first segment in the separated flow series and the wake width behind each separated flow 
segment is read into WKVV. The equivalent wake cross section area behind each segment is the wake 
width WKW multiplied by the projected length of the body segment axis on the tunnel z-axis. 


The Mode/ Wing. The wing lies in the z=ZWlNC! plane in the model coordinate system and is described 
by input data at the number of wing stations specified by NWS. The YWS array gives the y-coordinate 
of each station. If the symmetry option is selected (1SYM-=1) the stations given should apply to the 
half wing on one side of the y— *0 plane. The far field perturbations due to wing thickness and lift are 
approximated by representing the chordwise distributions of thickness and lift at each wing station by the 
first four members of a multipole singularity series located at the half-chord point of the wing station. The 
x-locations of the half-chord points are read into the XCW array. 

The multipole coefficients required to complete the wing input data at each wing station may be 
evaluated as follows. Let X l = x- XCW, let the thickness distribution t(x) and the lift distribution q(a:) 
be defined from the section leading edge Aj = —c/2 to the trailing edge X± = c/2. Further, express the 
thickness gradient as a = dt/dx. Then the thickness multipole coefficients are given by 

/•c/2 

QSO =- / adXi ~ t'pj? 

J-c/2 
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Note that the integrals of the form / CY[ l tLYj are the thickness coefficients of the series used by Rizk 
and Smithmeyer (ref. 5) which is applicable only for zero trailing-edge thickness. The present series 
can be evaluated from an “equivalent inviscid” thickness distribution in which t-gg represents the wake 
displacement thickness giving rise to wake blockage. 

The lift multipole coefficients at each wing station for a particular wing lift coefficient are given by 
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which are identical to those used in ref. 5. The vorticity in the above integrals is the lifting vorticity 
component at a particular wing lift coefficient normalized by the reference velocity. 

The Model Tail. The tail input quantities are evaluated in a manner completely analogous to that 
described above for the wing. Note, however, that the vorticity scaling is arbitrary for the tail because of 
the pitching moment matching procedure discussed below. 

Aerodynamic Coefficient Matching. The RANCOR program is able to make scaling and other adjustments 
to the test model data entered in record types 19 through 29 to be appropriate for the values of CLIFT, 
CD RAG and CMOM entered in record type 37. For matching lift coefficient, two sets of the wing lift 
multipole coefficients, corresponding to two different wing lift coefficients are entered in record type 27. 
The program evaluates the product C'jr • SREF for each set as two times the trapezoid rule integration 
over the full wing span of QC01 or QC02. A set of QGx appropriate for the input value of CLIFT is then 
calculated by linear interpolation (or extrapolation). The tail lift is accounted for in this process. 

The pitching moment match accounts for contributions of the wing, tail, and the line doublet segments 
used in the body representation. If a tail is not present, the input value of CMOM is matched by adding a 
constant to all values of QG1, thereby shifting the wing lift effects upstream or downstream as required. If 
a tail is present, the tail effectiveness is evaluated by appropriate integrations of QGOT and QG1T across 
the tail span and then all of the QGxT are scaled as required to achieve a match. If CMOM is input as 
99, all QGxT are scaled to zero and no adjustmant to the wing QG1 is made. 

A fairly crude approximation is used in matching the input value of CDRAG. It is assumed that the 
drag of the model representation is equivalent to complete loss of momentum in a wake having a cross 
section area equal to the net source strength of the model representation normalized by reference velocity. 
This simple approximation is believed to be acceptable for use in the PANCOIt program because with 
measured pressures used as wall boundary conditions, a change in net model source strength causes a 
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compensating change in wall flux with only minor effects on the assessed tunnel interference. The model 
net source strength is evaluated from the wing and tail trailing edge thicknesses and the point sources used 
in the body representation. After accounting for induced drag, any required adjustment is accomplished 
by a change in wing trailing edge thickness which is distributed across the span in proportion to the values 
of QS 1 . 

Sting Representation 

The input data form for the sting is analogous to that described for the model body in this Appendix. 
The sting station coordinates XSS and ZSS are expressed in the sting coordinate system having its origin 
at the sting center of rotation and oriented at a pitch angle TIIETS relative to the tunnel coordinate 
system. 

The sting is represented as a segmented inclined slender body by use of point sources and line doublets 
as described for the model body. In the case of the sting, however, the first point source, which would 
represent the growth in cross section area from zero to that of the first sting segment, is omitted. With 
this arrangement, a sting which continues the body lines behind the blunt base of a model body may be 
described with the first sting station located at the last body station. The combined representation is then 
equivalent to a sting fully replacing the wake behind the blunt body. If the sting is immersed in the wake 
behind a larger diameter model body base, the first sting station should be located farther downstream 
where the sting can be expected to start influencing the flow outside of the body wake. If the nose of the 
sting is exposed to the unshielded tunnel flow, the sting should be described with a dummy zero- volume 
segment placed ahead of the actual sting nose. 

It should be noted that the location of the model-sting interface is significant in that the interference 
perturbation at any point is calculated as the perturbation summed over all singularities except those 
included in the model representation. 
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Appendix D 
Sample Case 

The following data file is set up to perform an interference assessment of a Mach number 0.7 test 
point in the National Transonic Facility (NTF) at NASA Langley Research Center. The test section is 
about 8.2 feet square with solid side walls and six slots each in the top and bottom walls extending from 
tunnel station 0.5 to station 25.0. Pressure coefficients obtained from longitudinal rows of orifices halfway 
between slots on the top and bottom walls and a single pressure at station 23.125 on the side wall center 
line are used in the assessment problem. The test model is a generic subsonic transport model referred to 
as Pathfinder I. The volume of a vertical tail is accounted for by adjustments to the values of QSOT and 
QS1T at the horizontal tail center line. 

Fourteen panel networks are used to simulate the tunnel. Networks 1 and 14 provide the upstream 
and downstream closures respectively. Networks 2 and 9 are the doublet networks for the slotted top and 
bottom walls, networks 3 and 10 model the discrete wall slots, networks 4 and 11 are the IST=6 networks 
used to model the reentry flap region and wall step, and networks 5 and 12 are specified source strength 
networks used to represent wall divergence. Network 0 is the doublet network for the tunnel side wall and 
networks 7 and 8 are IST=2 prescribed source networks to represent the effective sidewall shape. The 
relatively sharp sidewall bend at the diffuser entrance is simulated by network 8 while network 7 provides 
an opportunity to represent boundary layer growth in the test section. The two phenomena are represented 
in separate networks to avoid unwanted upstream propagation of source gradients from the sharp bend 
into the test section. Network 13 approximates the sting support sector as a wedge-nosed plate. 

Sample Case Input Data File 

The record type number is shown in brackets in the left hand column to assist the user and is not a 
part of the input data file. 


[ii i. 

[ 2 ] 1110 

[3] 1 14 13.000 0. 

[4] 220010004 

5 25 0003 0044 12 444444 

5 15 20900 39 1 

452060306 
232010200 

6 25 0003 00444444444 

276020100 
266020100 

5 25 000 3 004444 12 4444 

5 15 90900 39 1 

459060306 
239010200 
220010100 
2200 -1 1 004 

[5] 100 l.E-10 1. 

[£] 0 . 0 . 0 . 0 . 0 . 0 . 0 . 0 . 

0 . 0 . 0 . 0 . 0 . 
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[?.] 

0. 

0. 0. 


vv 

0. 

0. 0. 

[8] 

3 

2 1 -10. 

[9] 

- A . 101042 A . 101042 

I ioj 

4.101042 

0. 

[8] 

2 

1 3 -4.101042 

[9] 

4.101042 

3.417505 2.050521 

[10] 

-10. 

-2. .5 



10.6 

12. 13. 



21.25 

22.5 23.75 

1 


30. 


[8] 

2 

1 3 -4.101042 

[9] 

4.101042 

3.417505 2.050521 

[10] 

.5 

2.0 3.6 

It 

13. 

14. 15. 

[8] 

2 

1 3 -4.101042 

(9 J 

4.101042 

2.597326 1.230313 

[10] 

20. 

21.25 22.5 

[8] 

2 

1 3 -4.101042 

[9] 

4.101042 

0. 

[10] 

.5 

25. 30. 

[8] 

3 

1 2 4.101042 

[9] 

4.101042 

2.460625 .820208 

[10] 

-10. 

-2. .5 



10.6 

12. 13. 



21.25 

22.5 23.75 



30. 



0 . 

0 . 

0 . 

0 . 

0 . 

0 . 



.683507 

0 . 



2.0 

3.6 

5.4 

7.2 

14. 

15. 

16.4 

18.2 

25. 

25.64 

26.28 

26.92 

.683507 

0 . 



5.4 

7.2 

9. 

10.6 

16.4 

18.2 

20. 

25. 

0 . 

23.75 

25. 




-.820208 

-2.460625 

-4.101042 


2.0 

3.6 

5.4 

7.2 

14. 

15. 

16.4 

18.2 

25. 

25.64 

26.28 

26.9 


[8] 

3 

1 2 4.101042 

[9] 

4.101042 

-4.101042 

[10] 

.5 

3.6 7.2 

[8] 

3 

1 2 4.101042 

[9] 

4.101042 

-4.101042 

[10] 

25.000 

25.640 26.280 

[8] 

2 

1 3 4.101042 

[9] 

0 . 

.683507 2.050521 

[10] 

-10. 

-2. .5 



10.6 

12. 13. 



21.25 

22.5 23.75 



30. 


[8] 

2 

1 3 4.101042 

[9] 

0 . 

.683507 2.050521 

[10] 

.5 

2.0 3.6 

It 

13. 

14. 15. 

[8] 

2 

1 3 4.101042 

[9] 

0 . 

1.230313 2.597326 

[10] 

20. 

21.25 22.5 

[8] 

2 

1 3 4.101042 

[9] 

0 . 

4.101042 

[10] 

.5 

25. 30. 

[8] 

3 

1 2 0. 

19] 

4.101042 

-4.101042 

[10] 

26.486 

29.642 

[8] 

3 

2 1 30. 

[9] 

4.101042 

-4.101042 

[10] 

4.101042 

0 . 


10.6 

13. 

15. 

30. 

26.920 

27.560 

30.000 


3.417505 

4.101042 



2.0 

3.6 

5.4 

7.2 

14. 

15. 

16.4 

18.2 

25. 

25.64 

26.28 

26.92 

3.417505 

4.101042 



5.4 

7.2 

9. 

10.6 

16.4 

18.2 

20. 

25. 

4.101042 

23.75 

25. 




0 . 


9. 

20 . 

27.56 


12 . 


9. 

20 . 

27.56 


9. 

20 . 

27.56 


12 . 


24 


.0214 

07276 

-.0691 

-.0987 

-.0862 

238 

1 23.125 

4.101042 

0 . 

1 

1 2.8 

2.734 

- 4.101042 

2 

1 4.5 

2.734 

- 4.101042 

3 

1 6.3 

2.734 

- 4.101042 

4 

1 8.1 

2.734 

- 4.101042 

5 

1 9.8 

2.734 

- 4.101042 

6 

1 11.3 

2.734 

- 4.101042 

7 

1 12.5 

2.734 

- 4.101042 

8 

1 13.5 

2.734 

- 4.101042 

9 

1 14.5 

2.734 

- 4.101042 

10 

1 15.7 

2.734 

- 4.101042 

11 

1 17.3 

2.734 

- 4.101042 

12 

1 19.1 

2.734 

- 4.101042 

13 

1 20.625 

2.734 

- 4.101042 

14 

1 2.8 

1.367 

- 4.101042 

15 

1 4.5 

1.367 

- 4.101042 

16 

1 6.3 

1.367 

- 4.101042 

17 

1 8.1 

1.367 

- 4.101042 

18 

1 9.8 

1.367 

- 4.101042 

19 

1 11.3 

1.367 

- 4.101042 

20 

1 12.5 

1.367 

- 4.101042 

21 

1 13.5 

1.367 

- 4.101042 

22 

1 14.5 

1.367 

- 4.101042 

23 

1 15.7 

1.367 

- 4.101042 

24 

1 17.3 

1.367 

- 4.101042 

25 

1 19.1 

1.367 

- 4.101042 

26 

1 20.625 

1.367 

- 4.101042 

27 

1 2.8 

0 . 

- 4.101042 

28 

1 4.5 

0 . 

- 4.101042 

29 

1 6.3 

0 . 

- 4.101042 

30 

1 8.1 

0 . 

- 4.101042 

31 

1 9.8 

0 . 

- 4.101042 

32 

1 11.3 

0 . 

- 4.101042 

33 

1 12.5 

0 . 

- 4.101042 

34 

1 13.5 

0 . 

- 4.101042 

35 

1 14.5 

0 . 

- 4.101042 

36 

1 15.7 

0 . 

- 4.101042 

37 

1 17.3 

0 . 

- 4.101042 

38 

1 19.1 

0 . 

- 4.101042 

39 

1 20.625 

0 . 

- 4.101042 

1 

1 2.8 

0 . 

4.101042 

2 

1 4.5 

0 . 

4.101042 

3 

1 6.3 

0 . 

4.101042 

4 

1 8.1 

0 . 

4.101042 

5 

1 9.8 

0 . 

4.101042 

6 

1 11.3 

0 . 

4.101042 

7 

1 12.5 

0 . 

4.101042 

8 

1 13.5 

0 . 

4.101042 

9 

1 14.5 

0 . 

4.101042 

10 

1 15.7 

0 . 

4.101042 

11 

1 17.3 

0 . 

4.101042 

12 

1 19.1 

0 . 

4.101042 

13 

1 20.625 

0 . 

4.101042 

14 

1 2.8 

1.367 

4.101042 



[ 15 ] 

15 

1 4.5 

1.367 

4.101042 







16 

1 6.3 

1.367 

4.101042 







17 

1 8.1 

1.367 

4.101042 







18 

1 9.8 

1.367 

4.101042 







19 

1 11.3 

1.367 

4.101042 







20 

1 12.5 

1.367 

4.101042 







21 

1 13.5 

1.367 

4.101042 







22 

1 14.5 

1.367 

4.101042 







23 

1 15.7 

1.367 

4.101042 







24 

1 17.3 

1.367 

4.101042 







25 

1 19.1 

1.367 

4.101042 







26 

1 20.625 

1.367 

4.101042 







27 

1 2.8 

2.734 

4.101042 







28 

1 4.5 

2.734 

4.101042 







29 

1 6.3 

2.734 

4.101042 







30 

1 8.1 

2.734 

4.101042 







31 

1 9.8 

2.734 

4.101042 







32 

1 11.3 

2.734 

4.101042 







33 

1 12.5 

2.734 

4.101042 







34 

1 13.5 

2.734 

4.101042 







35 

1 14.5 

2.734 

4. 101042 







36 

1 15.7 

2.734 

4.101042 







37 

1 17.3 

2.734 

4.101042 







38 

1 19.1 

2.734 

4.101042 






r 

39 

1 20.625 

2.734 

4.101042 





[ 16 ] 

10 

0 0 







[181 

0 

55 -.05 

4.101 

0 . 

26.95 

4.101 

0 . 




0 

55 0 . 

2.0505 

0 . 

27 , 

2.0505 

0 . 




1 

11 .4531 

0 . 

-.1771 

1.76555 

2.5 

.0408 




1 

11 2.3065 

0 . 

0 . 

2.922 

1 . 

.1763 




1 

23 - 2.17 

0 . 

-.1771 

3.33 

0 . 

-.1771 




1 

23 - 2.17 

.5 

-.1335 

3.33 

.5 

-.1335 




1 

23 - 2.17 

1 . 

-.0899 

3.33 

1 . 

-.0899 




1 

23 - 2.17 

1.5 

- . 0463 

3.33 

1.5 

-.0463 




1 

23 - 2.17 

2 . 

-.0027 

3.33 

2 . 

-.0027 


1 


1 

23 - 2.17 

2.5 

.0408 

3.33 

2.5 

.0408 


[ 19 ] 

22 

9 5 

0 0 






[ 20 ] 

.8301 

0 . 

0 . 

-.1 

.05 




[ 21 ] 

.8432 

0 . 

1.988 

.4783 

4.244 




[ 22 ] 

- 1.4633 

- 1.200 

-.9000 

-.6000 

-.3000 

0 . 

.1817 

.2292 

1 


.2767 

.3242 

.4287 

.5427 

.67577 

.7897 

.9037 

1.0177 



1.1317 

1.4 

1.75 

2.1 

2.4 

2.7033 



[ 23 ] 

0 . 

0 . 

0 . 

0 . 

0 . 

0 . 

-.0111 

-.01694 



-.01973 

-.02194 

-.02412 

-.02504 

-.02444 

-.02194 

-.01797 

-.01512 



-.01487 

0 . 

0 . 

0 . 

0 . 

0 . 



[241 

.01194 

.03954 

.05298 

.05412 

.05412 

.03258 

.00899 

.00920 



.00932 

.02072 

.02277 

.02660 

.02262 

.02228 

.02188 

.02170 



.04863 

.06039 

.05401 

.03914 

.03135 




[ 26 ] 

0 . 

.6567 

.95 

0 . 

0 . 

0 . 

0 . 




.24 

.6567 

.95 

0 . 

0 . 

0 . 

0 . 




.24 

.6567 

.95 

.005947 

-.08107 

-.06197 

-.04642 




.533 

.7717 

.7332 

.00387 

-.0427 

-.0252 

-.0144 




.8267 

.8871 

.5158 

.00355 

-.02167 

-.00940 

-.00396 




1.203 

1.08335 

.4463 

.00293 

-.0144 

-.00530 

-.00184 




1.5792 

1.27955 

.3768 

.002326 

-.01007 

-.00311 

-.00094 




1.9125 

1.4563 

.3153 

.001784 

-.00692 

-.00180 

-.00045 




2.2117 

1.6142 

.26 

.001318 

-.00470 

-.00102 

-.00022 
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[ 27 ] 

.108 

.00293 

.0143 

.00030 

.320 

-.0448 

.0369 

-.00727 



.105 

.00556 

.0115 

.00051 

.300 

-.0308 

.0277 

-.00466 



.105 

.00556 

.0115 

.00051 

.300 

-.0308 

.0277 

-.00466 



.097 

.00806 

.0067 

.00052 

.272 

-.0185 

.0165 

-.00183 



.085 

.00658 

.00258 

.00022 

.241 

-.0120 

.0074 

-.00070 



.070 

.00572 

.00134 

.00014 

.211 

-.0086 

. 0044 

-.00035 



. 055 

.00500 

.00073 

.00009 

.177 

-.0060 

.0027 

-.00018 



.040 

.00420 

.00035 

.00006 

.140 

-. 0038 ’ 

.0016 

-.00009 



.027 

.00320 

.00015 

.00003 

.102 

-.0017 

.0006 

-.00002 

128J 

0. 

2.412 

.5371 

.01981 

-. 14461 

-.00857 

-.00376 




.2068 

2.412 

.5371 

.002632 

-.01896 

-.00857 

-.00376 




.4101 

2.516 

.4261 

.002088 

-.01193 

-.00428 

-.00149 




.6133 

2.620 

.3152 

.001544 

-.00653 

-.00173 

-.00045 




. 8166 

2.723 

.2042 

.001001 

-.00274 

-.00047 

-.00008 


[ 29 ] 

.250 

-.075 

.0064 

-.0034 







.149 

-.025 

.0041 

-.00085 







.077 

-.0045 

.0013 

-.00012 







.049 

-.0072 

.00045 

-.00010 







.029 

-.0040 

.00020 

0. 





[ 30 ] 

34 

0 0 







[ 31 ] 

2.6162 

2.7670 

2.9266 

3.0953 

3.2737 

3.4625 

3.6621 

3.8733 



4.0966 

4.3327 

4.5825 

4.8468 

5.1262 

5.4217 

5.7343 

6.0649 



6 . 4146 . 

6.7844 

7.1756 

7.4272 

7.6967 

7.9855 

8.2948 

8.7762 

1 


9 . 3504 

10.5 

11.237 

11.975 

12.749 

13.488 

14.539 

15.592 

1 


16.644 

19.810 







[ 32 ] 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 



0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

0. 

[ 33 ] 

.00980 

.01069 

.01165 

.01272 

.01392 

.01524 

.01672 

.01938 


| 

.02440 

.03141 

.04014 

.05090 

.06419 

.08054 

.10055 

.11554 



.12271 

.13034 

.09132 

.11229 

.13809 

.16974 

.31293 

.45895 


f 

1.23863 

.9946 

1.1134 

1.2985 

1.3084 

1.7969 

1.6768 

1.5526 



4.4998 








[ 35 ] 

PATHFINDER I IN NTF 

' - SAMPLE 

CASE 





[ 36 ] 

21 . 

18 . 

399 . 

.79885 





[ 37 ] 

.79947 

2.58044 

.55950 

.03840 

-.07500 




[ 12 ] 

.05668 

.22867 

.47877 

.80678 





tt 

.05668 

.22867 

.47877 

.80678 





[ 38 ] 

.00186 

-.03484 







• 1 

.00195 

-.03487 







[ 39 ] 

.02742 








[ 40 ] 

.01448 

.01705 

.02099 

.01892 

.01554 

.01513 

.01179 

.01424 



.02047 

.02002 

.01145 

.00020 

-.00446 

.01616 

.01767 

.01724 



.01930 

.01646 

.01350 

.01437 

.01567 

.01575 

.01606 

.01096 



-.00295 

-.01267 

.01819 

.01940 

.02034 

.02034 

.01969 

.01916 



.01919 

.01939 

.01923 

.01773 

.01200 

-.00019 

-.01397 

.01938 



.01962 

.02062 

.01966 

.01495 

.00730 

.00103 

-.00130 

-.00047 



.00227 

.00308 

-.00312 

-.01235 

.01198 

.01847 

.02238 

.02406 



.01825 

.00486 

.00127 

-.00381 

-.00227 

.00818 

.00154 

-.01177 



-.02015 

.01887 

.01849 

.01968 

.02150 

.01646 

.00253 

-.00067 



-.00012 

.00214 

.00427 

.00168 

-.01100 

-.02303 
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Figure 1, Network superposition and boundary conditions to represent a slotted tunnel wall. 
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a. Standard output, MOVE(l)=l. 

Figure 3. Output printed by PANCOR program^ 
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. Continued. 
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a. Continued. 
Figure 3. Continued. 
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a. Continued. 
Figure 3. Continued. 
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a. Continued. 
Figure 3. Continued. 
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a. Concluded. 
Figure 3. Continued. 
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. Additional geometry output, MOVE(l)=2 or 4. 
Figure 3. Continued. 
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b. Continued. 
Figure 3. Continued. 
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Figure 3. Continued. 
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b. Concluded. 
Figure 3. Continued. 
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. Additional solution output, M0VE(l)=3 or 4. 
Figure 3. Continued. 
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Continued. 



PANEL CENTER OUTPUT NNE= 2 NL= 5 HL=25 
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